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ABSTRACT 

The meteorology of hot Jupiters has been characterized primarily with thermal measurements, but 
recent observations suggest the possibility of directly detecting the winds by observing the Doppler 
shift of spectral lines seen during transit. Motivated by these observations, we describe the types of 
Doppler signatures generated by the atmospheric circulation and show how Doppler measurements 
can place powerful constraints of the meteorology. We show that, depending on parameters, the 
atmospheric circulation — and Doppler signature — of hot Jupiters splits into two regimes. At weak- 
to- moderate stellar insolation, the day-night thermal forcing generates fast zonal jet streams from the 
interaction of atmospheric waves with the mean flow. In this regime, air along the terminator (as 
seen during transit) flows toward Earth in some regions and away from Earth in others, leading to a 
broadened or bimodal Doppler signature exhibiting distinct, superposed blue- and redshifted velocity 
peaks. Under more intense stellar insolation, however, the thermal forcing is so strong that it damps 
these planetary-scale waves, inhibiting their ability to generate jet streams. Strong frictional drag 
likewise damps these waves and inhibits jet formation. As a result, this second regime exhibits a 
circulation dominated primarily by high- altitude, day-to-night airflow along both terminators rather 
than longitudinally symmetric jets. This causes air to flow toward Earth along most of the terminator, 
leading to a predominantly blueshifted Doppler signature during transit. We present state-of-the-art 
three-dimensional circulation models including nongrey radiative transfer to quantify this regime shift 
and the resulting Doppler signatures; these models suggest that cool planets like GJ 436b lie in the 
first regime, HD 189733b is transitional, while planets hotter than HD 209458b lie in the second 
regime. Moreover, we show how the amplitude of the Doppler shifts place strong constraints on 
the strength of frictional drag in the upper atmospheres of hot Jupiters. If due to meteorology, the 
~2kms _1 blueshift inferred by Snellen and colleagues may require drag time constants as short as 
10 4 -10 5 seconds, possibly the result of Lorentz-force braking on this planet's hot dayside. 
Subject headings: planets and satellites: general, planets and satellites: individual: HD 209458b, 
methods: numerical, atmospheric effects 



1. INTRODUCTION 

To date, the exotic meteorology of hot Jupiters has 
been characterized primarily with thermal emissi on ob- 
servations, particularly infrared light curves (e.g. [K nut- 
son et al. 1120071 ICowan et al.| |2007| |Knutson et al.||2UU9| 
Crossfield et al.||2010D and secondary eclipse measure- 



ment s (e.g., Charbonneau et al. 2005 



20051). Together, these observations p' 



Deming et al. 



ace important 
constraints on the vertical temperature profiles, day- 
night temperature differences, and magnitude of day- 
night heat transport due to the atmospheric circulation. 
Moreover, in the case of HD 189733b and Ups And b, in- 
frared lightcurves indicate an eastward displaceme nt of 
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ture is a common out 



models, which generally exhibit fast eastward windflow 
at the e quator that displaces the thermal maxima to 
the east (Showman fc Guillot|2002j |Cooper fc Showman 



2005) IShowman eT"aT 



2009; [bobbs-ljixon Lin 
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|et al.||2010||Heng et al.||2011| |Showman fc Polvani|20TT 
In this way, the light curves provide information— albeit 
indirectly — on the atmospheric wind regime. 

Recent developments, however, open the possibility 
of direct observational measurement of the a tmospheric 
winds on hot Jupiters. |Snellen et al. (2010|) presented 



high-resolution groundbased, 2-fim spectra obtained dur- 
ing the transit of HD 209458b in front of its host star. 
From analysis of 56 spectral lines of carbon monoxide, 
they reported an overall blueshift of 2 ± 1 kms -1 relative 
to the expected planetary motion, which they interpreted 
as a signature of atmospheric winds flowing from dayside 
to night side towar d Earth along the pl anet's terminator. 
In a similar vein, |Hedelt et al.| ( |2011[ ) presented trans- 
mission spectra of Venus from its 2004 transit, in which 
they detected Doppler shifted spectral lines in the upper 
atmosphere, again seemingly the result of atmospheric 
winds. These observations pave the way for an entirely 
new approach to characterizing hot Jupiter meteorology. 

The possibility of characterizing hot Jupiter meteo- 
rology via Doppler provides a strong motivation for de- 
termining the types of Dopple r signatures generated by 



the atmo s pheric circulation. Miller- Ricci Kempton 



Rauscher (2012) took a first look at the ability of the 



atmospheric circulation to affect the transmission spec- 
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j3 effect) leads to the emergence of zonal j e ts, which of - 
ten dominate the circulation 



Fig. 1. — Schematic implications of atmospheric circulation 
for Doppler measurements of a hot Jupiter observed in transit. 
Planet's host star is depicted at right, planet (viewed looking down 
over north pole) is at the center, and Earth is at the left, (a): 
In the presence of zonal jets, air flows along latitude circles (col- 
ored arrows), leading to airflow toward Earth along one terminator 
(blue arrow) and away from Earth along the other (red arrow). A 
Doppler signature that is broadened, or in extreme cases may be 
bimodally split into blue- and redshifted components, results, (b): 
When zonal jets are damped, air flows primarily from day to night 
at low pressure, leading to airflow toward Earth along both ter- 
minators (blue arrows). A primarily blueshifted Doppler signature 
results. 

trum. Here, we continue this line of inquire to show how 
Doppler measurements can place powerful constraints 
on the meteorology. We show that the atmospheric 
circulation of hot Jupiters splits into two regimes — 
one with strong zonal jets and superposed eddies, and 
the other comprising predominant day-to-night flow at 
high-altitudes, with weaker jets — which exhibit distinct 
Doppler signatures. 

In Section |2j we present theoretical considerations 
demonstrating why two regimes should occur and the 
conditions for transition between them. In Section [3j 
we test these ideas with an idealized one-layer dynam- 
ical model. Section [4] presents state-of-the-art three- 
dimensional dynamical models of three planets — GJ 
436b, HD 189733b, and HD 209458b— that bracket a 
wide range of stellar irradiation and plausibly span 
the transition from jet to eddy-dominatecj^] at the low 
pressures sensed by Doppler measurements. Section [5] 
presents the expected Doppler signatures from these 
models, and Section [6] concludes. 

2. TWO REGIMES OF ATMOSPHERIC 
CIRCULATION: THEORY 

We expect the Doppler signature of the atmospheric 
circulation on hot Jupiters to fall into two regimes, illus- 
trated in Figure [T] 

2.1. Jet- dominated regime 

On rotating planets, the interaction of atmospheric 
turbulence with the anisotropy introduced by the merid- 
ional gradient of the Coriolis parameter (known as the 

4 Eddies refer to the deviation of the winds from their zonal 
average. 
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forcing and friction are weak, the heating of air parcels 
as they cross the dayside or nightside will be too small to 
induce significant day-night temperature variations; the 
dominant driver of the flow will then be the meridional 
(latitudinal) gradient in the zonal-mean radiative heat- 
ing. Such a flow will exhibit significant zonal symmetry 
in temperature and winds with the primary horizontal 
temperature variations occurring between the equator 
and the poles. For this regime to occur, the radiative 
timescale must be significantly longer than the timescale 
for air parcels to cross a hemisphere; the rotation rate 
must also be sufficiently fast, and the friction sufficiently 
weak. The speed and number of the zonal jets will de- 
pend on a zonal momentum balance between Coriolis 
accelerations acting on the mean-meridional circulation 
and eddy accelerations resulting from baroclinic and/or 
barotropic instabilities, if any. 

When the radiative forcing is sufficiently strong, as 
expected for typical hot Jupiters, large da y- night heat- 
ing con t rasts will occur. As shown by Showman & 
|Polvani (2011), such heating contrasts induce standing, 
planetary-scale Rossby and Kelvin waves. For typical 
hot Jupiter parameters, these waves cause an equator- 
ward flux of eddy angular momentum that drives a su- 



perrotating ( eastward) jet at the equator (Showman 
Polvani 2011[ ). This provides a theoretical explanation 



for the near-ubiquitous emergence of eastward equatorial 
jets in atmospheric circulation models of hot Jupiters. 

In these jet-dominated regime^] (Fig. [T^i), air along the 
terminator — as seen during transit — flows toward Earth 
in some regions and away from Earth in others. This 
leads to a Doppler signature where spectral lines are 
broadened, with minimal overall shift in the central wave- 
length. In extreme cases the Doppler signature may be 
split into distinct, superposed blue- and redshifted veloc- 
ity peaks. 

2.2. Suppression of jets by damping 

The presence of sufficiently strong radiative or fric- 
tional damping can suppress the formation of zonal jets, 
leading to a circulation that at high altitudes is dom- 
inated by day- to-night flow rather than jets that are 
quasi-symmetric in longitude (Fig. fib) . Here, we demon- 
strate the conditions under whicn the mechanisms of 



S howman fc Polvanirt20 1 lh are suppressed. 



Showm an & Polvani| Q2Q11| ) identified two specific 
mechanisms lor the emergence of equatorial superrota- 
tion in models of synchronously rotating hot Jupiters. 
We consider each in turn. 

2.2.1. Differential zonal wave propagation 



5 The interaction of eddies with the mean flow is generally re- 
sponsible for driving zonal jets, so eddies are almost never negli- 
gible to the dynamics, even when zonal jets are strong. Here, by 
"jet dominated" we do not mean that eddies are unimportant but 
rather simply that the resulting jets have velocity amplitudes that 
significantly exceed the amplitude of the eddies. 
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As described above, the day-night thermal forcing on 
a highly irradiated, synchronously rotating planet gener- 
ates standing, planetary-scale Rossby and Kelvin waves. 
The Kelvin waves straddle the equator while the Rossby 
waves exhibit pressure perturbations peaking in the mid- 
latitudes for typical hot Jupiter parameters. The (group) 
propagation of Kelvin waves is to the east while that 
of long Rossby waves is to the west; this differential 
zonal propagation induces an eastward phase shift of 
the standing wave pattern near the equator and a west- 
ward phase shift at high latitudes. The result is a pat- 
tern of eddy velocities (northwest-southeast in the north- 
ern hemisphere and southwest-northeast in the southern 
hemisphere) that causes an equatorward flux of eddy an- 
gular momentum. 

If the radiative or frictional timescales are significantly 
shorter than the time required for Kelvin and Rossby 
waves to propagate over a planetary radius, the waves are 
damped, inhibiting their zonal propagation and prevent- 
ing the latitude-dependent phase shift necessary for the 
meridional angular momentum fluxes. Therefore, this 
mechanism for generating zonal jets is suppressed when 
radiative or frictional damping timescales are sufficiently 
short. The Kelvin- wave dispersion relation in the primi- 
tive equations is 



Nk 



(1) 



where uo is wave frequency, k > and m are zonal 
and vertical wavenumbers, respectively, N is the Brunt- 
Vaisala frequency, and H is the scale height. The fastest 
propagation speeds occur in the limit of long vertical 
wavelength (m — » 0), which yields oj = 2NHk and thus 
phase and group propagation velocities of 2NH. The 
propagation time across a hemisphere is thus roughly 
a/NH, where a is the planetary radius. We thus expect 
this jet-driving mechanism to be inhibited when 



Trad < 



NH 



or r drag <C 



NH 



(2) 



10 8 m, H 



3 x 10 3 s 1 appropriate to a ver- 



For typical hot- Jupiter parameters (a 
400 km, and N 
tically isothermal temperature profile for a gravity of 
10 ms -2 and specific heat at constant pressure of 1.3 x 
10 4 Jkg^K- 1 ), we obtain a/NH ~ 10 5 sec. Thus, this 
mechanism should be inhibited when the radiative or 
drag time scales are much shorter than ^10 5 s. 

2.2.2. Multi-way force balance 

Even when the radiative timescale is extremely short 
and zonal propagation of Rossby and Kelvin waves is 
inhibited, an eddy-velocity pattern that promotes equa- 
torial superrotati on can occur under s o me c onditions. As 
pointed out by | Showman fc Polvani (2011) in the con- 
text of linear solutions, a three-way horizontal force bal- 
ance between pressure-gradient, Coriolis, and frictional 
drag forces can lead to eddy velocities tilted northwest- 
southeast in the northern hemisphere and southwest- 
northeast in the southern hemisphere if the drag and 
Coriolis forces are comparable. This occurs because drag 
generally points opposite to the wind direction, whereas 
the Coriolis force points to the right (left) of the wind 
in the northern (southern) hemisphere. When these two 



forces are comparable, balancing them with the pressure- 
gradient force requires that the horizontal wind rotates 
clockwise of the day-night pressure-gradient force in the 
northern hemisphere and counterclockwise of it in the 
southern hemisphere. In the limit of short r ra d when 
the horizontal pressure-gradient force points from day 
to night, these arguments imply that, at low latitudes, 
the eddy velocities tilt northwest-southeast in the north- 
ern hemisphere and southwest-northeast i n the southern 
hemisphere (see [Showman fc Polvani|2011| Appendix D). 

Even when frictional drag is too weak to play an im- 
portant role in the force balance, a similar three-way bal- 
ance between pressure-gradient, Coriolis, and advection 
forces can under appropriate conditions lead to veloc- 
ity tilts that promote equatorial superrotation. As air 
flows from day to night, the Coriolis force will deflect 
the trajectory of the airflow to the right of the pressure- 
gradient force in the northern hemisphere and to the left 
of it in the southern hemisphere. When the pressure- 
gradient force per mass, Coriolis force per mass, and ad- 
vective acceleration are all comparable, as expected un- 
der the Rossby number Ro ~ 1 conditions typical of hot 
Jupiters, then the deflection will be substantial. In the 
limit of short r ra d when the horizontal pressure-gradient 
force points from day to night, these arguments again 
imply that the eddy velocities tilt northwest-southeast 
in the northern hemisphere and southwest-northeast in 
the southern hemisphere. 

Now consider the effect of damping on this mecha- 
nism. To the degree that the radiative timescale is short 
enough for temperatures to be close to radiative equi- 
librium, radiative damping will not inhibit this mech- 
anism; however, strong frictional damping can prevent 
it from occurring. When the frictional force is much 
stronger than the Coriolis and advective forces, the hor- 
izontal force balance is no longer a mult i- way force bal- 
ance but rather becomes essentially a two-way balance 
between the pressure-gradient force and drag. In this 
case, winds simply flow down the pressure gradient from 
day to night. There is thus no overall tendency for pro- 
grade eddy- velocity tilts to develop, so the jet-pumping 
Reynolds stress, and the jets themselves, are weak. 

To quantify the amplitude of drag needed for this tran- 
sition to occur, consider a drag force per mass parame- 
terized by — v/rdrag? where v is horizontal velocity and 
Td m g is the drag time constant. The drag force dom- 
inates over the Coriolis force when Td ra g <C where 
f = 2Q sin <j) is the Coriolis parameter, Q is the planetary 
rotation rate (2tt over the rotation period), and <f> is lat- 
itude. Models of hot Jupiters predict flows whose domi- 
nant length scales are global, in which case the advective 
acceleration should scale as Z7 2 /a, where U is the char- 
acteristic horizontal wind speed. Drag will then domi- 
nate over the advection force when Td ra g <^ WIV^I) 1 / 2 , 
where | V4>| is the characteristic amplitude of the horizon- 
tal day-night pressure-gradient force on isobar^] given to 

6 The condition for dominance of drag over advection can be 
motivated as follows. When advection and drag are comparable, 
and both together balance the pressure-gradient force, it implies 
to order-of-magnitude that 
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order of magnitude by |V3>| ~ i?ATh or izA lnp/a, where 
is the specific gas constant, ATh or iz is the character- 
istic day-night temperature difference, and Amp is the 
range of \np over which this temperature difference ex- 
tends. For typical hot Jupiter parameters, both condi- 
tions imply drag dominance for Td rag <C 10 5 sec. When 
this condition is satisfied, the horizontal force balance 
is between the pressure-gradient and drag forces. As 
mentioned above, the resulting circulation at low pres- 
sure involves day-to-night flow with minimal zonal-mean 
eddy-momentum flux convergences in the meridional di- 
rection and weak zonal jets. 

2.2.3. Direct damping of jets by friction 

Frictional drag can also directly damp the zonal jets. 
A robust understanding of how drag influences the equi- 
librated jet speed — and hence a rigorous theoretical 
prediction of the amplitudes of drag needed to damp 
the jet — requires a detailed theory for the full, three- 
dimensional interactions of the global-scale planetary 
waves with the background flow, which is currently lack- 
ing. It is therefore not possible at present to provide 
a robust theoretical estimate of the amplitude of drag 
necessary to damp the zonal jets. Still, because the 
jets are fundamentally driven by global-scale waves that 
result from the day-night heating gradients ( jShowman] 
Polvani 2011), and because the radiative time con- 



stant increases rapidly with depth, we expect that the 
magnitude of zonal-mean acceleration of the zonal-mean 
zonal wind varies strongly with depth. These arguments 
heuristically suggest that the necessary frictional damp- 
ing times are less than a value ranging from 10 4 sec at 
low pressure (where the Doppler measurements sense) 
to 10 6 sec or more at high pressures, below the infrared 
photosphere. 

2.2.4. Recap 

When the jets and the waves that generate them are 
suppressed, the planet will tend to exhibit a large day- 
night temperature difference at low pressure, resulting 
in a large horizontal pressure gradient force between day 
to night that will drive a day-night flow (modified by 
the Coriolis effect) at low pressure. In this regime, air 
flows toward Earth along most of the terminator, leading 
to a predominantly blueshifted Doppler signature dur- 
ing transit. Mass continuity requires the existence of a 
return flow from night to day in the deep atmosphere 
(below the regions sensed by Doppler transit measure- 
ments). Because the density at depth is much larger 
than that aloft, the velocities of this return flow can be 
small. 

3. TEST OF THE TWO REGIMES WITH AN 
IDEALIZED MODEL 

We now demonstrate this regime transition in an ideal - 
ized dynamical model. As in |Showman fc Polvani] ( |2Q11[ ), 
we consider a two-layer model, with constant densities 
in each layer; the upper layer represents the stratified, 

These two relations yield Td rag ~ a/U and U ~ Td ra g|V$|, which 
together imply Td ra g ~ (a/| V^l) 1 / 2 . For drag time constants sig- 
nificantly shorter than this value, the drag force exceeds the ad- 
vection force. 



meteorologically active atmosphere and the lower layer 
represents the denser, quiescent deep interior. When the 
lower layer is taken to be infinitely deep and the lower- 
layer winds and pressure field are steady in time, the 
governing equations are the shallow-water equations for 
the flow in the upper layer: 



dv 

~dt 



■ gVh + /kxv = R- 



7"drag 



dh h eqK 
— + V • (vh) = — 

Ot T rad 



Q 



(4) 



(5) 



where v(A, (/>, t) is horizontal velocity, ft(A, 0, t) is the up- 
per layer thickness, A is longitude, t time, g is the (re- 
duced) gravity^ and d/dt = d/dt + v • V is the material 
(total) derivative. The term R in Eq. Q represents mo- 
mentum advection between the layers; it is — vQ/ft in 
regions of heatin g (Q > 0) and zero in reg ions of cool- 
ing (Q < 0). See Showman & Polvani (2011) for further 
discussion and interpretation ol the equations. 

In the context of a three-dimensional (3D) atmosphere, 
the boundary between the layers represents an atmo- 
spheric isentrope, and radiative heating/cooling, which 
transports mass between layers, is therefore represented 
as a mass source/sink, Q, in the upper-layer equations. 
We parameterize this as a Newtonian cooling that relaxes 
the thickness toward a radiative-equilibrium thickness 
ft, 

Here, we set 



eq (A, (/>), over a prescribed radiative time constant r ra d. 



ft eq (A,(/>) 



H 
H 



Aft eq cos A cos ( 



on the night side; 
on the dayside 

(6) 

where the substellar point is at (A, (f>) = (0°,0°). This 
expression incorporates the fact that, on the nightside, 
the radiative equilibrium temperature profile of a syn- 
chronously rotat ing hot Jupiter is constant (e.g., |Show-| 
|man et aL||2008| ), whereas on the dayside the radiative- 
equilibrium temperature increases from the terminator to 
the substellar point. An important property of Eq. |6| is 
that the zonal-mean radiative-equilbrium thickness, ft eq , 
is greater at the equator than the poles, reflecting the fact 
that a planet with zero obliquity (whether tidally locked 
or not) absorbs more sunlight at low latitudes than high 
latitudes. Note t hat Eq. fl6| differs from the formulation 



of ft eq adopted by Showman & Polvani (2011 ), where ft, 



was set to H + Aft eq cos A cos across the entire planet 
(dayside and nightside). 

In addition to radiation, we include frictional drag 
parameterized with Rayleigh friction, — v/rdrag, where 
Tdrag is a specified drag timescal e. The drag could resul t 
from vertical turbulen t mixing (|Li fc Goodlnan||2QTQ| , 
Lorentz-force braking (Perna et al. 2U10| ), or other pro- 
cesses. 

Our model formulation is identical to that described 



in Showman & Polvani (2011 Section 3.2) in all ways 
except tor the prescription ol ft e q(A, </>). 

7 The reduced gravity is the gravity times the fractional density 
difference between the two layers. For a hot Jupiter with a strongly 
stratified thermal profile, where entropy increases signficantly over 
a scale height, the reduced gravity is comparable to the actual 
gravity. 
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Parameters are chosen to be appropriate for hot 
Jupiters. We take gH = 4 x 10 6 m 2 sec -2 and set 
Ah eq /H = 1, implying that the radiative-equilibrium 
temperatures vary by order-unity from nightside to day- 
side. We also take Q = 3.2 x 10 _5 sec _1 and a = 
8.2 x 10 T m, implying a rotation period of 2.2 Earth days 
and radius of 1.15 Jupiter radii, similar to the values for 
HD 189733b. The radiative and frictional timescales are 
varied over a wide range to characterize the dynamical 
regime. 

We solved Equations (|4|-(5| in full spherical geome- 
try using the Spectral Transform Shallow Water Model 
(STSWM) of |Hack fc Jakob] ( [19921 ). The equations are in- 
tegrated using a spectral truncation of T 170, correspond- 
ing to a resolution of 0.7° in longitude and latitude (i.e., 
a global grid of 512 x 256 in longitude and latitude). All 
models were integrated until a steady state is reached. 

The solutions confirm our theoretical predictions of a 
regime transition. Figure [2] illustrates the equilibrated 
solutions for radiative time constants, T ra d, of 10, 1, 0.1, 
and 0.01 dayj jfor the case where drag is turned off (i.e., 
As expected, when the radiative time con- 
(10 days, panel (a)), jets dominate the circu- 
lation, with relatively weak eddies in comparison to the 
zonal-mean zonal winds. At intermediate values of the 
radiative time constant (1 and 0.1 days, panels (b) and 
(c)), the flow consists of strong jets and superposed ed- 
dies. At short values of the radiative time constant (0.01 
days, panel (d)), the jets are relatively weak — though 
not absent — and day-to-night eddy flow dominates the 
circulation. 

The dynamical behavior of this sequence can be un- 
derstood as follows. When the radiative time constant 
is long (Figure J2^a)) , the day-night thermal forcing is 
weak — air parcels experience only weak heating/cooling 
as they circulate from day to night — and the circula- 
tion is instead dominated by the equator-to-pole vari- 
ation in the zonal- mean h ec{ (i.e., by the zonal- mean ra- 
diative heating at low latitudes and cooling at high lati- 
tudes). At intermediate values of the radiative time con- 
stant (panels (b) and (c)), the day-night thermal forc- 
ing becomes sufficiently strong to generate a significant 
planetary wave response, and the eddy- momentum con- 
vergence induced by these waves generates e quatorial su- 
per rot at ion via the mechanisms identified by Showman & 
J 2011| ), particularly the differential zonal propa- 
gation of the standing Kelvin and Rossby waves. At short 
radiative time constant (panel (d)), such zonal propa- 
gation is inhibited but, as predicted by the theory in 
Section [2j the three-way force balance between pressure- 
gradient, Coriolis, and advection forces still generates 
prograde phase tilts in the velocities. Although, visually, 
the flow appears dominated primarily by day-to-night 
flow (Figure J2^d)) 9 these phase tilts still drive superrota- 
tion near the equator, and even at higher latitudes the 
zonal-mean zonal wind remains a significant fraction of 
the eddy wind amplitude. 
The transition from a regime dominated by jets to 

8 In this paper, 1 day is denned as 86400 sec. 

9 As described by Sho wman fc Polvani| ( |201l| , the coupling be- 
tween layers — specifically, mass, momentum, and energy exchange 
in the presence of heating/cooling — ensures that even cases with- 
out drag in the upper layer readily equilibrate to a steady state. 
All the models shown here are equilibrated. 
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Fig. 2. — Geopotential gh (orange scale) and winds (arrows) 
for the equilibrated (steady-state) solutions to the shallow-water 
equations (Equations [4j[5| in full spherical geometry assuming no 
upper-level drag (r^rag - ► oo) and T ra d=10, 1, 0.1, and 0.01 Earth 
days from top to bottom, respectively. Note the transition from 
a circulation dominated by zonally symmetric jets at long r ra d to 
one dominated by day-to-night flow at short r ra d. 



a regime dominated by day-to-night flow is even more 
striking when drag is included. Figure|3]shows a sequence 
of models with radiative time constants, T ra d, of 10, 1, 
0.1, and 0.01 days (as in Figure [2| but with Td ra g = 10r ra d 
in all cases. Overall, the trend resembles that in Figure[2j 
at long radiative time constants (panel (a)), the flow is 
dominated by high-latitude, highly zonal jets; at inter- 
mediate radiative time constants (panel (b) and (c)), the 
flow is transitional, exhibiting strong eddies associated 
with the standing planetary-scale wave response to the 
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day- night thermal forcing, and zonal jets driven by those 
eddies (cf Showman fc Polvani|2011 ); and at short radia- 
tive time constants (panel (d)), the Kelvin and Rossby 
waves are damped and the circulation consists almost 
entirely of day-to-night flow. As predicted by the the- 
ory in Section |5J drag in this case is strong enough to 
overwhelm the advection and Coriolis forces, leading to 
a two-way horizontal force balance between the pressure- 
gradient force and drag. As a result, there is no overall 
prograde phase tilt of the velocity pattern. The eddy 
forcing of the zonal-mean flow, and the jets themselves, 
are therefore weak. 

To better characterize the dominance of jets versus 
day-night flow, we performed integrations over a com- 
plete grid including all possible combinations of 0.01, 0.1, 
1, 10, and 100 days in r ra d and 0.1, 1, 10, 100, and oo 
days in Td rag . The typical amplitude of the jets can be 
characterized by the root-mean-square of the zonal-mean 
zonal wind variation in latitude: 



i r /2 

- / u(<t)) 2 d(t) 

K J -it/2 



1/2 



(7) 



where the overbar denotes a zonal average. To char- 
acterize the amplitude of the eddies, we adopt a metric 
representing the variation of the zonal wind in longitude: 
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and then determine the root-mean-square variations of 
this quantity in latitude: 
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The ratio of u rms to ^ e ddy,rms then provies a measure of 
the relative dominacne of jets versus day-night flow. 

These calculations demonstrate that jets dominate 
when both the radiative and drag time constant are long, 
and that day-to-night eddy flow dominates when either 
time constant is very short. This is shown in Figure [4j 
which presents the ratio w rms /w e dd y ,rms versus r ra d and 
Tdrag- The dashed curve in panel (b), corresponding to 
^rms/^eddy,rms = 1, demarcates the approximate transi- 
tion between regimes (jets dominate above and to the 
right of the curve, while day-night flow dominates below 
and to the left of the curve). Although extremely short 
values of either r ra d or Td ra g are sufficient to ensure eddy- 
dominated flow, the trend of the transition differs for the 
two time constants. When drag is weak or absent, r ra d 
must be extremely short — less than 0.1 day — to ensure 
eddy- rather than jet-dominated flow (Figure [4]). On the 
other hand, over a wide range of r ra d values, Td ra g need 
only be less than ~3 days to ensure eddy-dominated flow. 
The transition between jet and eddy-dominated regimes 
as a function of drag occurs more sharply when r ra d is 
large than when it is small. 

4. THREE-DIMENSIONAL MODELS 

We now demonstrate this regime shift in three- 
dimensional atmospheric circulation models including re- 
alistic radiative transfer. GJ 436b, HD 189733b, and HD 




-50 50 

longitude [deg] 



100 150 



Fig. 3. — Geopotential gh (orange scale) and winds (arrows) 
for the equilibrated (steady-state) solutions to the shallow-water 
equations (Equations |4j|5| in full spherical geometry assuming 
r drag — 10 r rad an d T ra d=10, 1, 0.1, and 0.01 days from top to 
bottom, respectively. Note the transition from a circulation domi- 
nated by zonally symmetric jets at long r ra d to one dominated by 
day-to-night flow at short r ra d. 

209458b are selected as examples that bracket a large 
range in stellar irradiation and yet are relatively easily 
observable, hence representing good targets for Doppler 
characterization. 

4.1. Model setup 

We solve the radiation hydrodynamics equations us- 
ing the Substellar and Planetary Atm ospheric Circula- 
tion a nd Radiation (SPARC) model of |Showman et al~ 
(2009). This model couples the dynamical core ol the 
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TABLE 1 
Model Parameters 



(days) 




100 



-1.193 



-0.521 0.152 0.824 

Log 10 (Zonal Jet/Eddy Ratio) 



1 .496 



Fig. 4. — Ratio w r ms/w e ddy,rmsj characterizing the ratio of jets 
to eddies, versus r ra d and Td ra g- Top: triple-dashed-dotted, 
dashed-dotted, dashed, dotted, and solid curves depict results 
from models with Td ra g values of 0.1, 1, 10, 100 days, and in- 
finite (i.e., no drag in the upper layer), respectively. Bottom: 
Two-dimensional representation of the same data. Colorscale 
depicts log 10 (w rms /ii e ddy,rms) with even log-spacing of the con- 
tour intervals. Values of u r ms/w e ddy,rms range from 0.06 (i.e., 
10 - 1.193) to 31 (j e ^ 10 1 - 496 ). The dashed black contour denotes 
Wrms/u e ddy,rms = 1 (i.e., log 10 (w rms /it eddy >rms ) = 0), giving an 
approximate demarcation between the jet- and eddy-dominated 
regimes. Jets dominate when both r ra d and Tdrag are long? an d 
eddies (essentially day-night flow) dominate when either time con- 
stant becomes very short. 



MITgcm (Adcroft et al. 12004), which solves the primi 



tive equations of meteorology in global, spherical geome- 
try, using pressure as a vertical coordinate, to the state- 
of-the-art, nongr ay radiative transfer scheme of |Marley| 
|fc McKayl ([1999), which solves the multi-stream radia- 
tive transfer equations using the correlated- k method to 
treat the wavelength dependence of the opacities. Here, 
we use a two-stream implementation of this model. To 
date, this is the only circulation model of hot Jupiters to 
include realistic radiative transfer, which is necessary for 
accurate determination of heating rates, temperatures, 
and flowfield. Gaseous opacities are calculated assum- 
ing local chemical equilibrium for a specified atmospheric 
metallicity, allowing for rainout of any condensates. We 
neglect any opacity due to clouds or hazes. 

Model parameters are summarized in Table [I] Al- 
though the atmo sphere of GJ 436b i s likely enriched in 



heavy elements (|Spiegel et al.| 2010 Lewis et al.||2010 



Madhusudhan & Seager 2011), solar metallicity repre- 



sents a reasonable baseline lor HD 189733b and HD 
209458b, and to establish the effect of differing stel- 
lar flux at constant metallicity we therefore adopt solar 





GJ 436b 


HD 189733b 


HD 209458b 


radius (m) 


2.69 x 10 7 


8.2396 x 10 7 


9.437 x 10 7 


rotation period (days) 


2.3285 


2.2 


3.5 


gravity (ms -2 ) 


12.79 


21.4 


9.36 


Pbase(bars) 


200 


200 


200 


Ptop(bars) 


2 x 10- 5 


2 x 10- 6 


2 x 10~ 6 


N r 


47 


53 


53 


metallicity (solar) 


1 and 50 


1 


1 



metallicity for the gas opacities in our nominal models of 
all three planets. To bracket the range of plausible metal- 
licities, we also explore a model of GJ 436b with 50 times 
solar metallicity. For HD 189733b and GJ 436b, we ne- 
glect opacity due to strong visible- wavelength absorbers 
such as TiO and VO, as TiO and VO are not expected 
for these cooler planets. For HD 209458b, secondary- 
eclipse measurements sug gest the presence of a strato- 
sphere (Knutson et al. 2008), and we therefore include 
TiO and VO opacity for this planet, which allows a ther- 
mal inversion due t o the large visible- wavelength opac- 
ity of these species ( |Hubeny et al. 2003| Fortney et al.| 
2008). Synchronous rotation is assumed for HD 189733b 
and HD 209458b (with a substellar longitude perpetu- 
ally at 0°); the rotation of GJ 436b, however, is assumed 
to be pseudqsychr onized with its slightly eccentric orbit 
( |Lewis et aL|2010| )F1 The obliquity of all models is zero 
so that the substellar point lies along the equator. 

Our nominal models do not include explicit frictional 
drag in the upper levels However, several frictional 
processes may be importa nt for hot Jupiters, including 
vertical turbulent mixing (Li & Goodman 2010), break- 
ing small-scale gravity waves ( Watkins &; Cho 2010), and 
magnetohydrodynamic drag flPerna et al.|2010p . The lat- 
ter may be particularly important for hot planets such as 
HD 209458b. Accordingly, we additionally present a se- 
quence of HD 209458b integrations that include frictional 
drag, which we crudely parameterize as a linear relax- 
ation of the zonal and meridional velocities toward zercF^l 
over a prescribed drag time constant, Td ra g- Within any 



10 This has only a modest effect on the results; synchronously 
rotating models of GJ 436b on circular orbits exhibit similar cir- 
culation patterns. 

11 All of the models include a Shapiro filter to maintain numer- 
ical stability. In some of the models, particularly those for HD 
209458b, we also include a drag term in the deep atmosphere be- 
low 10 bars; this allows the total kinetic energy of the model to 
equilibrate while minimally affecting the circulation in the upper 
atmosphere. Note that, for computationally feasible integration 
times (thousands of Earth days), models that include drag in the 
deep layers (but not at pressures less than ~10 bars) exhibit flow 
patterns and wind speeds in the observable atmosphere that are 
extremely similar to those in models that entirely lack large-scale 
drag. This is due to the fact that, even in such drag-free models, 
the wind speeds at pressures > 10 bars remain weak. For brevity, 
in this paper, we use the term "drag free" to refer to models lack- 
ing an explicit large-scale drag term, — v/Td ra g> m the observable 
atmosphere; nevertheless, it should be borne in mind that some of 
those models do contain drag in the bottommost model layers, and 
all of them include the Shapiro filter. 

12 In other words, we add a term — v/rd ra g to the horizontal mo- 
mentum equations, where v is the horizontal velocity. This simple 
scheme is called "Rayleigh drag" in the atmospheric dynamics lit- 
erature. 
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Fig. 5. — Global temperature (orange scale) and winds (arrows) 
for a sequence of 3D SPARC/MITgcm models at a pressure of 
0.1 mbar, where Doppler measurements are likely to sense. Pan- 
els (a), (b), (c), and (d) show solar-metallicity GJ 436b, 50 x so- 
lar GJ 436b, solar-metallicity HD 189733b, and solar-metallicity 
HD 209458b models, respectively. The vertical blue solid and red 
dashed lines show the location of the terminators 90° west and east 
of the substellar longitude, respectively. The substellar point is at 
longitude 0° in all panels. The models show a gradual transition 
from a circulation dominated by zonal jets (top) to one dominated 
by day-night flow at low pressure (bottom). 



given model, we treat Td rag as a constant everywhere 
within the domain. This is not a rigorous representa- 
tion of drag (for example, Lorentz forces will vary greatly 
from dayside to nightside and may act anisotropically on 
the zonal and meridional winds); still, the approach al- 
lows a straightforward evaluation of how drag of a given 
strength alters the circulation. 

For all three planets, we solve the equations on the 
cubed-sphere grid using a horizontal resolution of C32, 
corresponding to an approximate global resolution of 
128 x 64 in longitude and latitude. The lowermost N r — 1 



vertical levels are evenly spaced in log-pressure from an 
average basal pressure Pbase of 200 bars to a top pres- 
sure, Ptop? of 20/ibars for GJ 436b and 2 /ibars for HD 
189733b and HD 209458b. The uppermost model level 
extends from a pressure of p t0 p to zero. Our models of GJ 
436b and H P 18 9733b were originally presented in [Lewis 
et al. (2010) and Fortney et aL| ( 2010| ), respectively, while 
tor HL) 204958b we present new models here. These new 
integrations adopt 11 opacity bins in our correlated- A: 
scheme; detailed tests show that this 11-bin scheme pro- 
duces net radiative fluxes, heating rates, and atmospheric 
circ ulations very simila r to those of the 30-bin models 
(see Kataria et al.||2012 ). We integrate these models un- 
til they reach an essentially steady flow configuration at 
pressures < lbar. 

4.2. Results: nominal models 

Our nominal, three-dimensional models exhibit a fun- 
damental transition in the upper-atmospheric behavior — 
at pressures where Doppler measurements are likely to 
sense — as stellar insolation increases from modest (for 
GJ 436b) to intermediate (for HD 189733b) to large (for 
HD 209458b). This is illustrated in Figures [sj and |6j 
Figure [5] shows temperature and winds over the globe 
for models of GJ 436b, HD 189733b, and HD 209458b. 
Figure [6] presents histograms of the fraction of termi- 
nator arc length versus terminator zonal-wind speed for 
these same mo dels J^] which gives an approximate sense 
of how a discrete spectral line would be split, shifted, or 
smeared in frequency due to the Doppler shift of zonal 
winds along the terminator. To isolate the effect of dy- 
namics, the contribution of planetary rotation to the ve- 
locity is not included in Figure [6| although we will return 
to its effects subsequently. 

In the case of solar-metallicity GJ 436b (Figure [5^), 
the predominant dynamical feature is a broad superro- 
tating (eastward) jet that extends over all longitudes and 
in latitude almost from pole to pole. The jet exhibits sig- 
nificant wave activity, manifesting as small-scale fluctua- 
tions in temperature and zonal wind, particularly at the 
high latitudes of both hemispheres where the zonal- mean 
zonal winds peak. Nevertheless, the model exhibits little 
tendency toward a zonal-wavenumber-one pattern that 
would be associated with a predominant day-to-night 
flow. Save for small regions near the poles, the zonal 
winds at low pressure are everywhere eastward, imply- 
ing that, during transit, the zonal winds flow away from 
Earth along the leading limb and toward Earth along the 
trailing limb. This would lead to almost equal blueshifted 
and redshifted Doppler components, with a relative min- 
imum near zero Doppler shift (Figure [6h). 

Next consider HD 189733b (Figure^). The model 
again exhibits a superrotating equatorial jet, which is 
fast across most of the nightside — achieving eastward 
speeds of 4kms _1 — but slows down considerably over 
the dayside, reaching zero speed near the substellar 

13 For each model, we define a one-dimensional array Ui corre- 
sponding to the terminator velocity at 0.1 mbar versus terminator 
angle 0{ from to 360°. We define 20 velocity bins, equally spaced 
between the minimum and maximum terminator velocities from 
the array Uj. We then determine the fraction of the points in the 
Ui array that fall into each velocity bin. This is what is plotted in 
Figure [6] The qualitative results are similar when different choices 
are maae for the number of velocity bins. 
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Fig. 6. — Histograms showing the fraction of the full, 360° termi- 
nator, at a pressure of 0.1 mbar, moving at various speeds toward 
or away from Earth for the four models presented in Figure [5] 
Negatives values are toward Earth and positive values are away 
from Earth. Panels (a), (b), (c), and (d) show solar-metallicity 
GJ 436b, 50 x solar GJ 436b, solar-metallicity HD 189733b, and 
solar-metallicity HD 209458b models, respectively. This gives a 
crude sense of the Doppler splitting of a discrete spectral line. The 
models exhibit a transition from flows that exhibit both blue- and 
red-shifted components (top) to a flow whose Doppler signature 
would be predominantly blueshifted (bottom). 



point. Despite this variation, the zonal winds within 
the jet (latitudes equatorward of 60°) are eastward along 
both terminators. In contrast, the high-latitude zonal 
wind (poleward of 60° latitude) is westward along the 
western terminator and eastward along the eastern ter- 
minatoif 1 " 4 "! as expected for day-to-night flow. As seen 
during transit, along the trailing limb, the zonal winds 
flow toward Earth. Along the leading limb, they flow 
toward Earth poleward of 60° latitude and away from 

14 Eastern and western terminators refer here to the terminators 
90° of longitude east and west, respectively, of the substellar point. 



Earth equatorward of 60° latitude. Spectral lines would 
thus exhibit a broadened or bimodal character, with the 
blueshifted component considerably stronger than the 
redshifted component (Figure]^). HD 189733b is thus 
a transitional case between the two regimes discussed in 
Section [5J 

In the case of HD 209458b (Figure [5}l) , the strong su- 
perrotating jet continues to exist, but at and west of the 
western terminator it is confined substantially closer to 
the equator than in our GJ 436b or HD 189733b mod- 
els. Poleward of ~30° latitude on the western terminator, 
and everywhere along the eastern terminator, the airflow 
direction is from day to night. This implies that, as seen 
during transit, the trailing limb exhibits zonal winds to- 
ward Earth. The leading limb exhibits zonal winds that 
are toward Earth poleward of ~30° latitude and away 
from Earth equatorward of ^30° latitude. This would 
lead to Doppler shifts that are almost entirely blueshifted 
(Figure [6jl) . 

To summarize, these models exhibit a transition from 
a circulation dominated by zonal jets at modest insola- 
tion (GJ 436b) to one dominated by day-night flow at 
high insolation (HD 209458b). Qualitatively, this transi- 
tion matches well the predictions from our theory in Sec- 
tion 2 — as the stellar insolation increases, the effective 
radiative timescale decreases, and this damps the stand- 
ing planetary-scale Rossby and Kelvin waves, limiting 
their ability to drive a dominant zonal flow and lead- 
ing to a circulation comprised primarily of day-to-night 
flow at low pressure. We emphasize that the models in 
Figures [5}|6] do not contain frictional drag at the low pres- 
sures sensed by remote measurements, and so the only 
source of damping is the radiation (as well as the Shapiro 
filter, which exerts minimal effect at large scales). The 
models show that the regime transition occurs very grad- 
ually as stellar insolation is varied (Figures [5^-d). This 
is also consistent with theoretical expectations; as shown 
in Figure |4j when large-scale drag is absent, the radia- 
tive time constant must be decreased by over a factor 
of ~30 (from ~3 days to less than 0.1 day) to force the 
flow from the jet-dominated to eddy-dominated regime. 
Moreover, as discussed in Section |2j damping through 
radiation alone can inhibit differential zonal propagation 
of the planetary-scale waves, but the multi-way horizon- 
tal force balance between pressure-gradient, Coriolis, and 
advective forces can still produce prograde phase tilts 
near the equator. Thus, we still expect a narrow equato- 
rial jet over at least some longitudes. This can be seen in 
the nonlinear shallow- water solutions (see Figure [2^d)) 
and also explains the continuing existence of a narrow 
equatorial jet even for extreme radiative forcing in the 
3D models (Figure gd)). 

This regime transition manifests clearly in plots of ter- 
minator winds. Figure [7] shows the wind component pro- 
jected along the line of sight to Earth at the termina- 
tor for a sequence of models. Red represents velocities 
away from Earth (hence redshifted) while blue represents 
velocities toward Earth (hence blueshifted). The first, 
second, and third rows show our nominal models of GJ 
436b, HD 189733b, and HD 209458b, while the fourth 
row depicts our model of HD 209458b adopting frictional 
drag with Td ra g = 10 4 s. For GJ 436b, the leading limb 
is redshifted while the trailing limb is blueshifted. For 
HD 189733b, the redshifted portion — corresponding to 
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that only about ~10% of the high-altitude limb is red- 
shifted while -90% is blueshifted. In the HD 209458b 
model with frictional drag, the high-altitude winds are 
blueshifted over the entire terminator, completing the 
transition from a circulation dominated by jets to one 
dominated by high-altitude day-to-night flow. 

The regime transition discussed here is affected not 
only by the stellar insolation but also the atmospheric 
metallicity. Larger metallicities imply larger gaseous 
opacities (due to the increased abundance of H2O, CO, 
and CH4 ), and this moves the photosphere to low er pres- 
sures (cf |Spiegel et al ||2010| [Lewis et al.||2010| ), imply- 
ing that the bulk ol the starlight is then absorbed in a 
region with very little atmospheric mass. As a result, 
increasing the atmospheric metallicity enhances the day- 
side heating and nightside cooling per mass at the pho- 
tosphere even when the stellar insolation remains un- 
changed. The effects of this are illustrated in Figure [5}d, 
which shows a GJ 436b model identical to that i n Fig- 
ure [5fi exce pt that the metallicity is 50 times solar (Lewis 



ure[5kexce 
eTal!l[20T0| 



2010). Because of the greater absorption of stellar 



HD 209458b *w 

CO 




radiation at high levels, the atmosphere exhibits a large 
day-night temperature difference and significant zonal- 
wavenumber-one structure in the zonal wind, with strong 
longitudinal variations in the equatorial jet reminiscent 
of that in our HD 189733b model (compare Figures pfb) 
and (c)). Although eastward flow still dominates along 
most of the terminator, as in the solar-metallicity GJ 
436b model, the western terminator exhibits westward 
flow within ~30° latitude of the pole. Spectral lines as 
seen during transit still exhibit bimodel blue and red- 
shifts, but the blueshifts are now slightly more dominant 
(Figure ^p). 

Although we have focused on the existence of a regime 
transition in models with differing stellar fluxes, it is 
worth emphasizing that the same transition often occurs 
within a given model from low pressure to high pres- 
sure. Generally speaking, the radiative time const ants 
are short a t low pressure an d long at high pressure ( |Iro| 
|et al.||2005l|Showman et al.|2008[ ). The theory presented 
here therefore predicts that, as long as the incident stellar 
flux is sufficiently high and frictional drag is sufficiently 
weak, the air should transition from a day-to-night flow 
pattern at low pressure to a jet-dominated zonal flow 
at high pressures. Just such a patter n is seen in many 



Fig. 7. — Winds toward or away from Earth (colorscale, ms — 1 ) 
along the full, 360° terminator in a sequence of models as viewed 
during the center of transit. Colorscale is such that red (posi- 
tive) represents redshifted velocities while blue (negative) repre- 

sents blueshifted velocities. The radial coordinate represents log |Menou[|2010| |Heng et al.||201ip . Note, however, that 
pressure, and the plotted range is from 200 bars at the inside to 



published 3D hot Jupiter models (e.g. , |Cooper fc Show- 
|man|2005l|2006l|Showman et al.|2008, 2009; RauscherT 



2/ibar at the outside. The first, second, and third rows show our 
solar-metallicity nominal models of GJ 436b, HD 189733b, and 
HD 209458b, respectively. The fourth row shows our model of HD 
209458b where frictional drag is imposed with a drag time con- 
stant of 10 4 s. For each model, the left panel shows the winds 
alone, and the right panel shows the sum of the winds and the 
planet's rotation. From top to bottom, the transition from high- 
altitude velocities that have both blue and redshifted components 
to velocities that are entirely blueshifted is clearly evident. 



the equatorial jet — is confined to the low and midlati- 
tudes on the leading limb. As a result, at high altitudes, 
only about one-quarter of the limb is redshifted, while 
about three-quarters is blueshifted. For our nominal HD 
209458b model, the confinement of the equatorial jet to 
low latitudes on the leading limb is even stronger, such 



if the incident stellar tiux is sufficiently low (and the 
drag is very weak), the atmosphere may be in a regime 
of jet-dominated flow throughout; on the other hand, if 
frictional drag is sufficiently strong, jets may be unable 
to form at all, and the atmosphere may be in a regime 
of day-night flow aloft with a very weak return flow at 
depth. 

4.3. Results: influence of drag 

We now consider the effect of frictional drag in 3D 
models. As discussed in Section[2j sufficiently strong fric- 
tional drag (i) damps the standing planetary-scale waves 
that are the natural response to the day-night heating 
gradient, and (ii) drives the horizontal force balance into 
a two-way balance between pressure-gradient and drag 
forces, both of which inhibit the development of prograde 
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Fig. 8. — Temperature and winds at 0.1 mbar pressure in models of HD 209458b with frictional drag. Left: a model with a drag time 
constant of 3 x 10 4 s. Right: a model with a drag time constant of 10 4 s. Top panels show temperature in K (orange scale) and winds 
(arrows). Bottom panels show zonal wind in ms — 1 . In the case with weaker drag, an equatorial jet extends partway across the nightside, 
but the jet is damped out in the case with stronger drag. 



phase tilts in the eddy velocities and in turn the pumping 
of zonal jets, and, finally (iii) directly damps the zonal 
jets. Thus, we expect that an atmosphere with suffi- 
ciently strong frictional drag will lack zonal jets and that 
its circulation will instead consist primarily of day-to- 
night flow at high altitude, with return flows at depth. 
Figure [8] illustrates this for solar-metallicity models of 
HD 209458b where Rayleigh drag is implemented with 
time constants of 3 x 10 4 s (left column) and 10 4 s (right 
column). As predicted, the air at 0.1 mbar flows directly 
from day side to nightside over both terminators. The 
model with Td ra g = 3 x 10 4 s exhibits a remnant equa- 
torial jet on the nightside that extends from the eastern 
terminator to the antistellar point. Because the stronger 
frictional drag damps it out, the model with Td ra g = 10 4 s 
lacks such a jet, and the flow exhibits only modest asym- 
metry (due to the /3 effect) between the western and 
eastern terminators. Doppler lines would be entirely 
blueshifted in both cases (Figure^). 

Friction affects not only the qualitative circulation 
regime (e.g., existence or lack of zonal jets) but also the 
speed of the high- altitude flow between day and night. 
Figure [9] shows the root-mean-square zonal wind speeds 
at the terminator for a sequence of HD 209458b mod- 
els with differing drag time constants. All of these runs 
are in the same regime as the model in Figure [5J where 
day-to-night flow dominates at low pressure. When Td ra g 
is sufficiently long, the flow speeds are independent of 
the drag time constant, but they start to decrease when 
Tdrag is sufficiently short. When drag is absent in the up- 
per atmosphere, our HD 209458b model equilibrates to 
an rms terminator wind speed of 5.2 kms -1 at 0.1 mbar, 
decreasing with depth to 3.8, 2.6, and 1.9 kms -1 at 1, 
10, and 100 mbar, respectively. As shown in Figure [9j 
the addition of weak drag (rd ra g = 3 x 10 5 s) exerts only 
a modest effect on the day-night flow speeds at pressures 
< 10 mbar. Drag time constants Tdrag 

< 10 5 s start to 
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Fig. 9. — Steady-state root-mean-square wind speeds at the ter- 
minator versus frictional drag time constant from a sequence of HD 
209458b models including drag. For each 3D model, performed for 
a given drag time constant, the root-mean-square wind speeds — 
calculated along the terminator — are shown at pressures of 0.1, 1, 
10, and 100 mbar. For these models the speeds generally represent 
day-night flow. The equilibrated speeds depend significantly on 
the drag time constant and, within a given model, on pressure. 



matter significantly in the upper atmosphere, however; 
for example, for Td ra g = 10 4 s, the rms terminator speeds 
are 2.1, 1.1, 0.6, and 0.2 kms -1 at 0.1, 1, 10, and 100 
mbar — significantly less than the equilibrated speeds in 
the absence of upper-level drag. 

The above results suggest that the amplitude of the 
observed Doppler shift can place constraints on the 
strength of f rictional drag in the upper atmospheres of 
hot Jupiters. jSnellen et al.| ( |2010| )'s inference of winds on 
HD 209458b is tentative but, at face value, suggests wind 
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speeds toward Earth of 2±1 kms -1 . Snellen et al. (2010) 
suggest that their measurements are sensing pressure lev- 
els of 0.01-0.1 mbar. At these levels, the winds in our 
models equilibrate to ~4-6kms _1 when drag is weak, 
and Figure [9] shows that reducing the wind speed to 
2 kms -1 requires drag time constants potentially as short 
as ~10 4 s. This hints that strong frictional drag processes 
may operate in the atmosphere of HD 209458b. But 
caution is warranted: Figure [9] also demonstrates that 
the rms terminator wind speeds also depend strongly 
on pressure within any given model; therefore, making 
robust inferences about drag amplitudes from observed 
Doppler shifts requires extremely careful and accurate 
estimates of the pressure levels being probed. This may 
be a challenge, at least until the composition and hence 
wavelength-dependen t opacity of hot Jup iters are better 
understood. If the |Snellen et al.| ( |2010| ) measurements 
are actually sensing deeper pressures of ^10 mbar, say, 
then explaining their 2-kms _1 signal would require little 
if any drag in the observable atmosphere. 

The qualitative dependence of terminator wind speed 
on the drag time constant — illustrated in Figure [9] — can 
be understood analytically. To order of magnitude, the 
horizontal pressure gradient force in pressure coordinates 
between day and night can be written i?ATh rizA lnp/a. 
This is balanced by some combination of advection, of 
magnitude U 2 /a, Coriolis force, of magnitude fU, and 
drag, of magnitude ?7/ r drag- Drag will dominate when 
Tdmg < / _1 and when /7/r drag > U 2 / 'a, which requires 
r drag <^ ( a /\ V^l) 1 / 2 , equivalent to the requirement that 
Tdrag <^ ^/(^AThoriz A Inp) 1 / 2 . As long as these con- 
ditions are satisfied, we can balance drag against the 
pressure-gradient force. Solving for Td ra g then implies 
that the amplitude of drag necessary to obtain a wind 
speed U is 

Ua 

Tdrag ~ i?AT horiz Alnp (10) 

Inserting parameters appropriate to the 0.1 mbar level 
on HD 209458b (a ~ 10 8 m, R ~ 3700 J kg" 1 K" 1 , 
AThoriz ~ 1000 K, and Amp ~ 5), and adopting U ~ 
2 kms -1 motivated by the Snellen et al. measurement of 
HD 209458b, we obtain Tdrag ~ 10 4 s. This value agrees 
well with the strength of drag needed in our 3D model 
intregrations to achieve a speed of 2 kms -1 at the 0.1 
mbar level (leftmost black triangle in Figure [9| . 

5. TRANSMISSION SPECTRUM CALCULATIONS 

To quantify the implications for observations, in this 
section we present theoretical transmission spectra from 
our 3D models demonstrating the influence of Doppler 
shifts due to atmospheric winds. These spectra illus- 
trate how the dynamical regime shifts described in the 
preceding sections manifest in transit spectra. 

5.1. Methods 

We have previously developed a code to compute the 
transmission spectrum of transiting planet atmospheres, 
which we extend here to include Doppler shifts due to at- 
mospheric winds. The first generation of the code, which 
used one-dimensional atmosp heric pressure-tempera ture 
(p-T) profiles, is described in |Hubbard et al.| (|2001) and 



Fortney et al. 



2003). In Shabram et al. (2011) the one- 



dimensional (lb) version of the code was well- validated 
against the analytic transmission atmosp here model of 



iLecavelier Pes Etangs et al. (2008). In Fortney et al. 
(2010J) we implemented a method to calculate the trans 
mission spectrum of fully 3D models. 

The calculation of the absorption of light passing 
through the planet's atmosphere is based on a sim- 
ple physical picture. One can imagine a straight path 
through the planet's atmosphere, parallel to the star- 
planet-observer axis, at an impact parameter r from this 
axis. The gaseous optical depth tq, starting at the ter- 
minator and moving outward in one direction along this 
path, can be calculated via the equation: 



r' cr(r')n(r') 



2^/2 



dr f , 



(11) 



where r' is the distance between the local location in the 
atmosphere and the planetary center, n is the local num- 
ber density of molecules in the atmosphere, and a is the 
wavelength-dependent cross-section per molecule. Later 
we will discuss the role of winds leading to a Doppler 
shifted o~ away from rest wavelengths. We assume hydro- 
static equilibrium with a gravitational acceleration that 
falls off with the inverse of the distance squared. The 
base radius is taken at a pressure of 10 bars, where the 
atmosphere is opaque, and this radius level is adjusted 
to yield the best fit to observations, where applicable. 
Here we define the wavelength-dependent transit radius 
as the radius where the total slant optical d e pth re aches 
0.56, following Lecavelier Des Etangs et al. (2008). Ad- 
ditional deta il and description can be found in |Fortney 



et al. ( |2010j ), as the 3D setup here is the same as de- 
scribed in that paper. 

For any particular column of atmosphere, hydrostatic 
equilibrium is assumed, and we use the given local p-T 
profile to interpolate in a pre-tabulated chemical equi- 
librium and opacity grid that extends o ut to 1 /ibar. 
The equilibrium chemistry mi xing ratios (Lodders 1999 



Lodders & Fegley 2002, 2006) are pa ired with the opac 

1_l |2008| ) to generate pres- 
ength- dependent absorp- 



ity database (|Freedman et al. 
sure - , temperature-, and wave, 
tion cross-sections that are used for that particular col- 
umn. For a different column of atmosphere, with a dif- 
ferent p-T profile, local pressures and temperatures will 
yield different mixing ratios and wavelength-dependent 
cross-sections. 

We include the Doppler shifts due to the local atmo- 
spheric winds and planetary rotation when evaluating the 
opacity at any given region of the 3D grid. At high spec- 
tral resolutio n, rotation tends to cause a broadening of 
spectral lines ( S piegel et al.|2007| ) , while the atmospheric 
wind speeds lead to absorption featu res that are Doppler 
shifted from their rest wavelen gths ( Snelle n" et al.|[20To 
|Miller-Ricci Kempton fc Rauscher|2012| ). The cross sec- 
tion a is not evaluated at the rest wavelength, Ao, but 
rather at the Doppler shifted value, A, found via 



A = A (l 



(12) 



where v\ os is the line-of-sight velocity — including both 
rotation an d atmospheric winds — and c is the speed of 
light. The |Snellen et al.| ( |2010[ ) observations were per- 
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formed at a resolving power of R ~ 10 5 . For additional 
clarity in presentation, we have computed opacities and 
transmission spectra at R = 10 6 . In practice we inter- 
polate within our R = 10 6 opacity database to yield the 
correct a for every height in the atmosphere, on every 
column, given the calculated velocities at every location. 
This is done at 128 locations around the terminator, and 
angles of ± ~ 20 degrees onto the day and night sides 
of the terminator. Note that for simplicity we do not 
include the Doppler shift due to orbital motion, and we 
are therefore essentially evaluating the transmission spec- 
trum at the center of the transit for a planet with zero 
orbital ecc entricity. The effect of orbital motion was con- 
sidered by Miller- Ricci Kempton & Rauscher (2012). 



5.2. HD 209458b Cases and the Role of Drag 

We now turn to a detailed analysis of the HP 2 09458b 
models in the vicinity of the |Snellen et al. (2010) obser- 
vations of HD 209458b near 2.2 /im. The Doppler shift 
was not measured across all wavelengths, but only within 
the narrow CO lines, since flat transmission spectra (cor- 
responding to the continuum between the spectral lines) 
yield no leverage on the Doppler shifts. T hese peaks are 
all of nearly the same strength (see, e.g. |Snellen et al.| 
2010, supplemental online material), so they probe very 
similar heights in the atmosphere. In Figure [10| we have 
computed the transmitted spectrum at R = 10 6 , 1300 
wavelengths, from 2.3080 to 2.3011 /im, for three mod- 
els of HD 209458b: a model with no drag in the upper 
atmosphere (top left panel), a weak-drag model with a 
drag time constant of 3 x 10 5 s (middle left panel), and a 
strong-drag model with a drag time constant of 1 x 10 4 s 
(bottom left panel). 

For each of the three drag cases we calculated the 
model planet radius vs. wavelength for four variants of 
the dynamical model. One use rest-wavelength values 
(black, this corresponds to a reference case where winds 
are assumed to be zero), one includes atmospheric dy- 
namics only (magenta, ignoring rotation but using the 
full 3D winds), one includes only rotation (green, ig- 
noring dynamics), and in orange is the full model, in- 
cluding both dynamics and rotation. The transmission 
spectra in Figure [l0| are somewhat difficult to interpret. 
Therefore we have also calculated the cross-correlation, 
compared to the rest wavelength model, across the 1300 
wavelengths in our calculation. 

There are several aspects of note for these plots. Start- 
ing in the upper right, the drag free HD 209458b model, 
we see that the self cross-correlation is strongly peaked 
at m s -1 , as expected. The dynamics-only model 
shows winds that peaked at —2500 ms -1 (meaning a 
blueshift), with strong winds (generally at the equator) 
reaching beyond —5000 ms -1 . As seen in the dynamics 
output, there is also a component of redshift winds, tak- 
ing up a relatively small fraction of the terminator, which 
range from ~0 to 5000 ms -1 . Interestingly, for this drag- 
free case, the cross-correlation curve of the rotation-only 
model is not symmetric about the zero- velocity point. 
This asymmetry only appears in models with a strong 
leading/trailing hemispheric temperature contrast. It 
appears to be due the trailing (hotter) hemisphere hav- 
ing a larger scale height, and therefore more prominent 
absorption features. The full model, including dynamics 



and rotation, has a broader peak than the dynamics- 
only model due to rotational broadening. The peak 
—5000 ms -1 velocities from dynamics and —2000 ms -1 , 
from rotation, lead to velocities on the trailing hemi- 
sphere's equator of —7000ms -1 . The full model is not 
merely just a broadened dynamics-only model due to the 
asymmetric rotational component. 

The weak-drag case (middle right) has a more con- 
strained atmospheric flow, which is generally day-to- 
night, with a much reduced super-rotating jet. Veloc- 
ities from dynamics are peaked more narrowly around 
-2200 ms -1 (red curve). The rotational component is 
close to symmetric, due to the small leading/ trailing tem- 
perature difference. The full model looks much like the 
dynamics-only model, broadened due to planetary rota- 
tion. 

The strong-drag case (lower right) has a very con- 
strained circulation, with a relatively uniform day-to- 
night flow around the entire planet, with little sign 
of an equatorial jet. The dynamics velocities are 
strongly peaked at —1000 ms -1 and the very small lead- 
ing/trailing temperature contrast lead to a symmetric ro- 
tational component. The full model looks very much like 
a broadened version of the dynamics-only model. The 
slightly higher peak on the right side of the full model 
is due to a slight asymmetry in the velocities around 
— 1000 ms -1 in the dynamics-only model. 



Overall, a comparison of the models in Figure 10 high- 
lights the possibility of inferring the amplitude ofdrag in 
the atmosphere of HD 209458b. The no-drag, weak-drag, 
and strong-drag models (Figure 10) exhibit peaks in the 
cross-correlation centered at —4, —2.5, and — lkms -1 , 
respectively. We emphasize that these values are the 
quantitative result of our rigorously calculated transmis- 
sion spectra from our fully coupled 3D models (and are 
not, for example, simply the velocity at some assumed 
pressure of the 3D models). At face value, the results 
in Figure [To] indicate that models with a d rag time con- 
stant o f 10^10 6 s provide a better fit to the Sne llen et al.| 
(2010) observations than models with no drag in the up- 



per atmosphere. Th is is consistent with the inferences 
drawn in Section [431 

The possibility of drag in the atmosphere of HD 
209458b is particularly interesting in light of recent sug- 
gestions that thermal ionization of alkali metals at high 
temperature can lead t o Lorentz forces that act to brak e 
the atmospheric winds ( Perna et al.||2010{ |Menou||2012[ ). 
In the regime of day-to-night tiow, air at the termma- 
tor has just crossed much of the dayside, and so the 
wind speed at the terminator is predominantly deter- 
mined by drag on the dayside rather than the nightside. 
Secondary-eclipse observations indicate that HD 209458b 



phere with temperatur es poten- 
2000 K dKnutson et al.| ( |2008[ ); s ee also 



exhibits a dayside stratos 
tially reaching 

our Figur e |5| . Based on the scaling relations in Perna 
et al. ( 2010[ )7such high temperatures should lead to very 
short drag times, potentially consistent with the infer- 
ences on drag drawn here. 

5.3. Comparing Three Different Planets 

Figure [TT] allows us to diagnose the different atmo- 
spheric dynamics and Doppler shift signatures of HD 
209458b, HD 189733b, and GJ 436b. All cases are drag- 
free. The top row is the same HD 209458b model de- 
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Fig. 10. — Left panels: High spectral resolution (R ~ 10 6 ) planet radius versus wavelength for three models of HD 209458b with various 
drag strengths. Solid black is the "Rest" model calculated without a Doppler shift. Magenta includes Doppler shifts due to atmospheric 
dynamics only. Green includes Doppler shifts due to planetary rotation only. In orange, the "Full Model" includes both dynamics and 
rotation. Right panels: Cross-correlations for the planetary transmission spectra shown at left. As drag becomes stronger one generally 
finds slower and more peaked wind speeds. See text for additional details. 



scribed in the top row of Figure [lOj The velocities for 
HD 189733 (middle panels) from dynamics are relatively 
similar to those of HD 209458b, but lacking a very high 
velocity component. The planet's rotation period of 2.2 
days is only 63% of the period of HD 209458b, meaning 
HD 189733b has a significantly larger rotational veloc- 
ity component. The asymmetry shown is also due to 
a relatively large leading/trailing temperature contrast. 
The full model has a smaller, high velocity peak at -5000 
m s _1 due to the strong blue-shifted peak of rotational 
velocity. 

The GJ 436b spectrum clearly shows a transition to 
a different regime of atmospheric dynamics. The flow 
is dominated by a wide super-rotating jet, with little 
flow purely being day to night. This manifests itself in 
the dynamics-only model as being somewhat symmetric, 
with a slightly higher cross-correlation peak on the blue- 
shifted side. The rotational component is nearly sym- 
metric, owning to relative temperature homogenization 



of the planet. The magnitude of the rotational velocities 
are quite small because the planet has a small radius. 
The full model shows little Doppler shift. The spectrum 
plot at left shows very little difference between all of the 
models, other than that they have been Doppler broad- 
ened compared to the rest model. 

6. DISCUSSION 

Recent observations suggest that atmospheric winds 
on hot Jupiters can be directly inferred vi a the Doppler 
shift of spectral lines seen during transit (Snellen et al. 
2010 ). Motivated by these observations, we have shown 
that the atmospheric circulation of hot Jupiters divides 
into two regimes depending on the strength of the radia- 
tive forcing and frictional drag, with implications for the 
Doppler signature: 

• Under moderate stellar fluxes and weak to mod- 
erate drag, atmospheric waves generated by the 
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Fig. 11. — Left panels: High spectral resolution (R ~ 10 6 ) planet radius versus wavelength for drag- free models of HD 209458b, HD 
189733b, and GJ 436b. Solid black is the "Rest" model calculated without a Doppler shift. Magenta includes Doppler shifts due to 
atmospheric dynamics only. Green includes Doppler shifts due to planetary rotation only. In orange, the "Full Model" includes both 
dynamics and rotation. Right panels: Cross-correlations for the planetary transmission spectra shown at left. As the incident flux becomes 
smaller, from HD 209458b, to HD 189733b, to GJ 436b, peak wind speeds become smaller, and the planet's dynamics become more 
dominated by a super-rotating jet. See text for additional details. 



day-night thermal forcing interact with the mean 
flow to produce fast east-west (zonal) jets. In this 
regime, air along the terminator flows toward Earth 
in some regions and away from Earth in others, 
leading to blue-shifted and red-shifted contribu- 
tions to the Doppler signature seen during transit. 
Depending on the speed of the winds relative to 
the planetary rotation, as well as the variation of 
zonal winds in latitude and height, this will cause 
Doppler lines observed during transit to be broad- 
ened or, in extreme cases, split into distinct, super- 
posed blue-shifted and red-shifted velocity peaks. 

• Under extreme stellar fluxes and/or strong fric- 
tional drag, however, the radiative and/or fric- 
tional damping is so strong that it damps these 
waves and inhibits jet formation. At the low 
pressures sensed by transit measurements, the at- 



mospheric circulation then involves a day-to-night 
flow, with a return flow at deeper levels. In this 
regime, the airflow at levels sensed by transit mea- 
surements is toward Earth along most or all of the 
terminator, leading to a predominantly blue-shifted 
Doppler signature of spectral lines observed during 
transit. 

We presented a theory predicting this regime transition, 
and we confirmed its existence and explored its prop- 
erties in one-layer shallow-water models and in three- 
dimensional models coupling the dynamics to realistic 
non-grey radiative transfer. We then presented detailed 
radiative transfer calculations of the transit spectra ex- 
pected from our 3D mode l s in th e 2- /am spectral region 
observed by Snellen et al. (2010); these calculations can 
help to guide future observational efforts. 

We also showed that, in the second regime described 
above, the speed of the day-night windflow depends on 
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the amplitude of the drag at the low pressures sensed 
by transit measurements. Under relatively weak drag, 
the wind speeds at the terminator of our HD 209458b 
models reach ~4-6kms _1 depending on altitude and 
forcing conditions. Under strong drag, the wind speeds 
are slower. Interestingly, at the low pressures sensed by 
transit observations, the drag must be relatively strong — 
with effective drag time constants of ~10 6 s or less — to 
reduce the speeds by a significant fraction. Our mod- 
els of HD 209458b without significant drag in the upper 
atmosphere produce peak cross-correlations of the tran- 
sit spectrum corresponding to blueshifts of ~3-7kms _1 ; 
this exceed s, albeit marginally, t he ~2kms _1 blueshift 



inferred by 
our models 



Snellen et al. 
that agree 



fl2010| ). 

test with the 



On t he other hand, 
Snellen et al.l 



(2010) inference — where the peak cross-correlations of 
the transit spectrum lie at ~l-3kms _1 — exhibit drag 
timescales of 10 4 -10 6 s. This suggests, tentatively, that 
frictional drag may be important on the dayside of HD 
209458b. An attractive possibility is that the dayside 
of HD 209458b is sufficiently hot for partial ionization 
to occur, leading to Lorentz- fo rce braking of the winds 
( Perna et al.|2 010; M enou|2012| ). Regardless, these mod- 
els demonstrate that, m principle, measurements of the 
Doppler shift of spectral lines can place constraints on the 
amplitude of drag in the atmospheres of hot J up iters. 

In light of this issue, we note that our results diffe r 
from those of Miller- Ricci Kempton fc Rauscher] ( |2012[ ) , 
whose obtained peak cross-correlations in the transmis- 
sion spectra corresponding to blueshifted velocities of 
about — 2.5kms _1 and — lkms -1 in models without and 
with drag, respectively. Their velocity shifts for their 
drag-free case are significantly slower than the veloc- 
ity shifts we obtain of — 4kms _1 for our drag-free HD 
209458b model. A significant difference in the two stud - 



ies is that, in Miller- Ricci Kempton & Rauscher (2012), 



the heating/cooling in the thermodynamic energy equa- 
tion was determined using a simplified Newtonian relax- 
ation schem e based on that presented in |Cooper fc Show^] 
man ( |2005[ ); in contrast, our dynamical models are fully 
coupled to non-grey radiative transfer, from which the 
radiative heating/cooling rates are calculated. This may 
lead to a quantitative difference in the radiative heating 
rates and hence three-dimensional wind structure. Fu- 
ture work may shed light on the discrepancies between 
the results. 

It is worth mentioning that, within each of the two 
broad dynamical regimes studied in this paper, there 
may lie additional subregimes involving important tran- 
sitions between dynamical mechanisms. For example, in 
the regime of zonal jets, we have emphasized the develop- 
ment of equatorial superrotation by standing, planetary- 
scale Rossby and Kelvin waves induced by th e day- night 
thermal forcing (jShowman fc Polvani |20lT ). However, 
when the stellar flux is lower than considered in most hot- 
Jupiter models, the importance of the day-night thermal 
forcing decreases and the equator-to-pole heating gradi- 



ent becomes dominant. Baroclinic instabilities can then 
occur, particularly when the rotation rate is fast, and 
these may lead to multiple mid-latitude east-west jets, 
with an equatorial jet that may be of either sign de- 
pending on the details. A transition analog ous to this is 
evident in models of GJ 436b presented by Lewis et al. 
(2010). We will explore such dynamical transitions fur- 
ther in future work. 



It is also wort h discussing the proposal of |Mon- 
talto et al. (2011), who suggested that t he ~2 kms -1 
blueshift inferred by |Snellen et al.| ( |2010[ ) results not 
from atmospheric winds but from planetary orbital mo- 
tion due to an eccentric orbit. The radial velocity 
at the time the planet crosses the line of sight to 
Earth is RVo = Ke cosvu /\Jl — e 2 where e is the ec- 
centricity, w is the longitude of periastron, and K = 
|27rG/P] 1 / 3 M^sinz/(M^ + ra p ) 2 / 3 = 1.47 x lO^sec" 1 
is a constant, which we have evaluated for HD 209458 
parameters. Here, P is the orbital period, G is the grav- 
itational constant, i is the orbital inclination, and M* 
and m p are the mass of the star and planet, respectively. 
Explaining a 2-km se c -1 blueshift would th us require 

(120111) point out 



Montalto et al. 



that e cos zd « 0.014. 
that the eccentricity itself is rather poorly constrained; 
however, what matters is not eccentricity alone but the 
combination ecosw. A key p oint, apparently not ap- 
preciated by Montalto et al. (2011), is that e cos zu is 
tightly constrained by observations of transit and sec- 
ondary eclipse. Observations of t he relative timing of 
transit and secondary eclipse from Deming et al. (2005) 
show that ecoszu < 0.002 at 1-a. Observations from 



Knutson et al.| (|2008) and |Crossfield et al.| ( 2012[ ) place 
even tighter upper limits on ecosw; the latter study 
yields ecostu = 0.00004 ± 0.00033, corresponding to a 
3-<j upper limit of the orbit-induced Doppler shift of 
140ms _1 at the center of transit. This appears to rule 
ou t any orbital explana tion for the Doppler shift inferred 
b yJSnellen et al.| ( |2010| ). 

Finally, while we have emphasized the wind patterns 
and implications for transit Doppler measurements, the 
dynamics described here also predict a regime transition 
in the temperature structure that may be important in 
explaining thermal observations from light curves and 
secondary eclipses. The shallow- water models in Fig- 
ures [2] and [3] and the 3D models in Figure [5j show that 
the now tends to a state with small longitudinal temper- 
ature variations when radiation and friction are weak, 
whereas the day-night temperature differences become 
large when either radiation or friction become strong. 
Our models therefore predict a transition from small to 
large fractional day-night temperature differences at the 
infrared photosphere as stellar flux increases from small 
to large. We will explore this issue further in future work. 



This research was supported by NASA Origins and 
Planetary Atmospheres grants to APS. 
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